Low-energy direct muon transfer from H to Ne^'^^, S^^+ and Ar 
using two-state close-coupling approximation to 
Faddeev-Hahn-type equation 

Renat A. Sultanov^'^* and Sadhan K. Adhikari^''' 
^Instituto de Fisica Teorica, Universidade Estadual Paulista, 
014-05-900 Sao Paulo, Sao Paulo, Brazil 
"^Department of Physics, Texas A & M University,^ 
College Station, Texas 77843 
(February 2, 2008) 



Abstract 



We perform three-body calculation of direct muon-transfer rates from 
thermalized muonic hydrogen isotopes to bare nuclei Ne^'^^, S^^^, and Ar^^"^ 
employing integro-differential Faddeev-Hahn-type equations in configuration 
space with two-state close-coupling approximation scheme. All Coulomb po- 
tentials including the strong final-state Coulomb repulsion are treated exactly. 
A long-range polarization potential is included in the elastic channel to take 
into account the high polarizability of the muonic hydrogen. The transfer 
rates so calculated are in good agreement with recent experiments. We find 
that the muon is captured predominantly in the n = 6, 9, and 10 states of 
muonic Ne^'^"'', S"*^^^, and Ar-^^^, respectively. 
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I. INTRODUCTION 



A theoretical description of three-body charge-transfer reactions in physics is both chal- 
lenging and interesting and is a subject of active research as the study of this funda- 
mental problem provides a testing ground of the underlying physical model. This is even 
more so when the mass of the exchanged charged particle is not negligible compared to 
the other two and when there is a strong final-state Coulomb interaction. The first feature 
invalidates the commonly-used simplifying Born-Oppenheimer-type approximation scheme 
and calls for a detailed three-body description. The second feature demands a proper 
dynamical treatment of the final-state Coulomb interaction. 

The direct muon (/i) transfer rates at low energies from the muonic-hydrogen isotopes 
^H^ (p/i) and ^H^ (rf/i) to bare nuclei X^"*" of large charge Z: 

(H^)i. + X^+ ^ (X^)(^-i)+ + H+ (1) 



have been the subject of many experimental investigations p|-[17|, where H stands for the 



hydrogen isotope proton (p) or deuteron (c?) and X-^^ stands for the target nuclei. Investi- 
gations of the low-energy muon-transfer reactions in collisions of muonic hydrogen with 
nuclei X^"*" are of importance for the muon- catalyzed fusion cycle |18|. Although, there 



have been many experiments of muon transfer with nuclei of large charge Z, the theoretical 
investigations are limited to nuclei YJ"^ with Z = 1 [0] and Z = 2 [§. The theoretical 
description of these reactions at low energies has both the complicating features mentioned 
above and becomes extremely complicated as the charge Z of the nuclei increases. This is 
due to the strong electromagnetic field of the nuclei simultaneously responsible for a large 
initial-state polarization and a strong final-state Coulomb interaction. Also, the large mass 
of the transferred muon compared to the electron leads to additional difficulties compared 
to the electron-transfer reaction where a Born-Oppenheimer-type approximation is efficient 
10]. In addition, a large number of open channels in the muon-transfer reactions even at 
zero energy complicates the theoretical treatment. It is difficult to incorporate these effects 
properly in a dynamical three-body calculation. This is why there are no published work 
to date on a three-body dynamical calculation of these muon-transfer rates for nuclei with 
charge Z > 3. These aspects demand a careful three-body quantum mechanical treatment 
with exact inclusion of the final-state Coulomb interaction. 

The recent theoretical activities in charge transfer are centered around problems with 
much weaker polarization and Coulomb interactions compared to those encountered in muon 
transfer involving bare nuclei like argon or sulphur. In electron- hydrogen-atom deuteron- 
hydrogen-atom 0, positron- hydrogen-atom 0, electron-positronium-atom ||^, hydrogen- 
positronium-atom and H-/i-H [0 systems the final-state Coulomb interactions are zero 
compared to 17e^/ p in the case of p-/x-Ar^^+ considered here, where e is the electronic charge 
and p is the radial separation in the final state. In proton-deuteron and H-/i-He systems 
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there is a small final-state Coulomb repulsion, e.g., e^/p. The initial-state polarization is 
also much smaller in those studies The large charge of the bare nuclei complicates 

substantially the analysis of the present muon transfer problems. 

Recently, we presented a theoretical formalism for the study of such muon-transfer reac- 
tions using two-component integro-differential Faddeev-Hahn-type equations in config- 
uration space with close-coupling approximation scheme |^ , and applied it to the study of 
muon transfer from muonic hydrogen isotopes to He^^ and Li^"*" , and C^"*" and O^^ . 
The transfer rates \z calculated there are in good agreement with experiments, whenever 
available. The formalism of Ref. seems to be appropriate for the study of the problem 
of muon transfer from muonic hydrogen isotopes to nuclei with large charge. 

Among such transfer reactions, muon transfers from muonic hydrogen isotopes to Ne^°^, 
S^^+, and Ar^^+ have been the center of active interest to different experimental groups. 
There have been several experiments for muon transfer from p/i to argon which yielded 
transfer rates centered around the following five different values: (16.3 ± 0.9) x 10^°s~^ 
@], (12.0 ± 1.9) X lO^^s-i 0, (14.6 ± 1.4) x lO^^s^^ [|n|], (35 ± 6) x lO^^s^^ and 



(98 ± 15) X 10^°s^^ [T^. For muon transfer from to argon the experimental transfer rates 
are (8.6 ± 0.4) x 10^°s~^ and 9.4 xlO^°s~^ lH]. The experimental situation in the case 
of argon is quite controversial with widely different values for the transfer rates, specially in 
the case of p/i. Despite this intense experimental activity in these muon-transfer reactions 
there are no quantum dynamical calculations valid at low energies for these transfer rates. 
In view of the above interest and controversy, in this paper we undertake the challenging 
three-body study of direct muon transfer from p/i and to argon using the formulation 
of Ref. P2[. To test our approach, in addition, we apply it to study muon transfer from 



hydrogen isotopes to neon and sulphur where experimental results are available. In the case 
of Sulphur the experimental transfer rate from p/i is 8.9 x 10^°s^^, whereas that from d^ is 
11 X 10^*^s^^ [|l^]. In the case of neon the experimental transfer rate from d^i is 10.1 x 10^°s~^ 



Although there is a very large number of open channels in these problems, for a given 
nuclei the muon is transferred predominantly to a few (muonic) atomic labels of the heavy 
nuclei X'^+ p3|-[25[]. For example, it was first noted by Gershtein and reconfirmed later 
1 2^ that the muon is captured mostly in the n = 4 states of C^^, and n = 5 states of 0^+ 



17| , |22[| . Also these transfers take place mostly to the final muonic-atomic states with low 
angular momenta and transfer rates are negligible for muonic atomic states with angular 
momenta I > 2. Semi classical description of these muon transfer reactions has been very 



useful in explaining many qualitative and quantitative features p^ , p6| . Using a semiclassical 
model based on potential curves of the two-centered p/iX^"*" system, Holzwarth et al. PB[ 
demonstrated that muon transfer to F^"*" takes place essentially to the n = 6 level of the 
(F^)^"*" system. They also showed that transfer to the nearby levels of the (F^)^"*" atom is 
negligible compared to the n = 6 level. From similar consideration of semiclassical barrier 
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penetration, Haff et al. derived the following formula for the state n of the (X^)"^^ 
atom to which the muon from is predominantly transferred in reaction (|l|) 

where In denotes the integral part. For C^'*', O^"*" and F^"'' this formula leads to n = 4, 
5, and 6, respectively, in agreement with calculation |p3| , |25|l . The formula is expected 
to work even better for heavier nuclei Ne^°^, S^^+ and Ar^^+ where the channelization to a 
final specific state n should be more perfect with muon transfer to nearby n states heavily 
suppressed. 

The correct dynamical formulation should include all open transfer channels and we in- 
cluded them in a previous study on muon transfer with light nuclear targets |21]. However, 
it is quite impossible now to treat even this reduced number of open transfer channels in a 
quantum calculation with heavier targets due to convergence difficulties in the presence of 
the large final-state Coulomb interaction mentioned above. Hence, in the present treatment 
we use a two-channel model to calculate transfer to a single final state, where we include 
the elastic and one transfer channel. Different sets of equations are used for the different 
final states. Eventually, the total transfer rate is calculated by summing the different con- 
tributions. After calculation we find that in both cases (p/i and d^) the muon is captured 



predominantly in the n = Q state of Ne 



lo^ 



n 



9 state of 8^^+, and 



n 



10 state of Ar^^+ 



in complete agreement with formula (0) and in this pioneering theoretical study we present 
results for these cases. The transfer is highly suppressed to other values of n and higher I 
states. 

In the Faddeev-Hahn-type equation |jT9|, the wave function is broken up into components 
with proper asymptotic behavior in different physical channels. Consequently, these wave- 
function components are much smoother functions of configuration-space variables than the 
Schrodinger wave function which is the sum of all these components. Hence, this approach 
simplifies the solution procedure and the correct asymptotic behavior of the solution in dif- 
ferent channels can be incorporated easily. In addition, these equations allow us to introduce 
explicitly a polarization potential in the initial channel. All Coulomb potentials including 
the strong final-state Coulomb repulsion are treated exactly (without approximation or pa- 
rameters) in this formalism. The effect of strong polarization of the muonic hydrogen by 
the bare nuclei is accounted for by a polarization potential with a cut off parameter. By a 
proper inclusion of the polarization potential in the intermediate region, the present transfer 
rates are found to be essentially independent of this parameter. The correct inclusion of the 
final-state Coulomb interaction has the advantage of building in the proper asymptotic be- 
havior of the wave function in a low-order close-coupling type approximation p7| , p8| . Hence 
as in Ref. p2| we make a two-state close-coupling approximation to the Faddeev-Hahn-type 
equation in the present study and find that a numerical solution using the present scheme 
leads to very encouraging agreement with recent experimental transfer rates. 
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The theoretical framework for the present study is based on the formahsm developed in 
Refs. p0|-[2^ which was used for the study of muon transfer from muonic hydrogen atoms 
to lighter charged nuclei H'*", He^"*", Li^+, C^"*", and O®^. In the dynamical equations in Refs. 



21,22] the final-state Coulomb interaction in the transfer channel is treated exactly without 



approximation. In addition, as in Ref. ||22|, here we explicitly include a polarization potential 
in the elastic channel. The presence of the strong Coulomb interaction and the associated 
large polarization make the present calculational scheme far more complicated numerically 
compared to those of Refs. . In a coupled- channel approach for atomic processes, the 



coupling to infinite number of p-wave states is responsible for generating the polarization 
potential p9|j30|] . As it is impossible to include all such states in a numerical scheme, the 
commonly accepted procedure is to replace these coupling terms by the polarization potential 
as in Ref. . 



In Sec. II we present a brief account of the theoretical formulation. In Sec. Ill we report 
the present numerical results and finally, in Sec. IV we present some concluding remarks. 



II. THEORETICAL FORMULATION 



A detailed account of the theory has already appeared p0|-|22|. For the sake of complete- 
ness we report here a summary of the theoretical development. Here we describe the dynam- 
ical equations we use based on the close-coupling approximation to Faddeev-Hahn-type two- 
component equations |^|. We use units e = h = = 1, where m^(e) is the muonic mass 



(charge), and denote the heavy nucleus X^~^ by 1, the hydrogen isotope (s) by 2 and muon by 
3. Below the three-body breakup threshold, following two-cluster asymptotic configurations 
are possible in the system 123: (23) — 1 and (13) — 2. These two configurations correspond 
to two distinct physical channels, denoted by 1 and 2, respectively. These configurations are 
determined by the Jacobi coordinates {rj3,Pk)- ^13 = "^3— "^i, P2 = (^3+"^i^i)/(l+"^i)— "^2, 
^23 = "^3 — ^2, pi = (r3 + m2r2)/(l + ^2) — ^i, where r^, rrii {i = 1, 2, 3,) are coordinates 
and masses of the particle i, respectively. 

Let us introduce the total three-body wave function as a sum of two components 

^(ri,f2,f3) = ^l(f23,Pi) + ^2(ri3,P2) (3) 

where ^1/1(^23, pi) is quadratically integrable over the variable r23, and "^2(^13, P2) over ri^. 
The components and \l/2 carry the asymptotic boundary condition for channels 1 and 2, 
respectively. The second component is responsible for pure Coulomb interaction in the final 
state. These components satisfy the following set of two coupled equations 

[E - {Ho + ^23(^23)) - f/pol(Pl)]^l(r23,Pi) = [(^3(^23) + Vl2(fi2)) - Uc (p2)]^2(ri3,P2) (4) 

[E - (Ho + V^i3(ri3)) - Uc (p2)]^2(ri3,P2) = [(V^i3(ri3) + ^12(^2)) - t/poi(pi)]^i(r23, Pi) (5) 
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where E is the center-of-mass energy, Hq the total kinetic energy operator, Vij{rij) the pair 
potential {i j = 1,2,3), Uq the final-state Coulomb interaction given by 

f/c(p.) = (6) 

with Z the charge of the heavy nuclei and Z'{= 1) the charge of the hydrogen isotope. Here 
f/poi is the polarization potential given by p5| 

and zero otherwise. The value of the cut-off parameter A has to be chosen appropriately. By 
adding equations @ and @ we find that they are equivalent to the Schrodinger equation. 

Distortion potentials are useful in phenomenological models of scattering in atomic 
physics Although, unnecessary in a complete solution of the Schrodinger equation. 



f/poi(pi) = for pi > A (7) 



they facilitate the numerical effort in a simplified model. We have included the proper po- 
larization potential in the initial channel. Although, a polarization potential exists in the 
final rearrangement channel, the most important interaction in this channel is the Coulomb 
repulsion which has been exactly included in our description. 

Because of the strong final-state Coulomb repulsion, it is very difficult to solve the 
multichannel model equations based on Eqs. (@) and (||). Hence, for solving @) and (H) we 
expand the wave function components in terms of bound states in initial and final channels, 
and project these equations on these bound states. The expansion of the wave function is 
given by 

^i(r23,p'i) ^ ^^^R[fli\r23\)/^n (8) 
Pi 

^2(ri3,p'2) ^ ^^^RZM^D {YciP2) ® nf,s)}oo (9) 

P2 

where nlC are quantum numbers of the three-body final-state, fii = 77131712/ {m^ + 7712), 
P2 = 77137711/(7713 + mi), Yim's the spherical harmonics, R^ni^f^-{\r\) the radial part of the 
hydrogen-like bound-state wave function for reduced mass /ij and charge Z, fil\pi) and 
fnic{P2) the unknown expansion coefficients. This prescription is similar to that adopted 
in the close-coupling approximation. After a proper angular momentum projection, the set 
of two-coupled integro-differential equations for the unknown expansion functions can be 
written as 



fil\pi) = gi^/{2jC + l) dp2fSkp2) 

JO \ \r23\ \ri2\ ) 

xpiP2C°Lri„(z/2,7r)/v^ (10) 



6 



ikl^^r + ^ - ^^^^ - 2M,Ucip2) fliP2)=92^i2C + l) 

I dpifil\pi) f dusinuRlf {\fis\) ( , , , , . , , 
"'0 Jo \ |ri3| |ri2| / 

xR[fXi\r23\)p2PiC'^oioYiUi^i,n)/VA^ . (11) 



z z 



Here k^^^ = ^2Mi{E - 4^) = ^2M2{E - E^n^) with Mf ^ = m^^ + (1 + ma)"^ and 

= m^^ + {l+mi)^^ , E^'' is the binding energy of pair (j3) and gj = AnMj/'y^ (j = 1, 2), 
7 = 1 — 17111712/ {{I + nil) (1 + 7712)), Ccoim the Clebsch-Gordon coefficient, L the total angular 
momentum, u the angle between the Jacobi coordinates pi and p2, i^i the angle between r23 
and pi and 1^2 the angle between ris and p2. 

To find unique solution to ( p!0D and ([TTl) , appropriate boundary conditions are to be con- 
sidered. We impose the usual condition of regularity at the origin f[l^ (0) = and f^c{0) = 0. 
Also for the present scattering problem with 1 + (23) as the initial state, in the asymptotic 
region, two solutions to (jT^) and (|TT|) satisfy the following boundary conditions 

fil\pi) ~ sin(fci^)pi)+i^nCOs(M'Vi) (12) 



P2-*+00 

where K^'' are appropriate coefficients. For scattering with 2 + (13) as the initial state, we 



fSkp2) ~ vW^^i"2Cos(^!'V2-r//2A:fhn2A;fV2-7r/:/2) (13) 

no — ^4-00 * 

have the following conditions 

fil\pi) - vW^^2"/cos(tfVi) (14) 
fSkP2) ~ sin(M'V2-r//2A;f)ln2A;fV2-vr£/2) 

+K^^ cos(fcf V2 - v/2k^i^ In 2A;f ^2 - vr/:/2) (15) 

where Vi [i = 1, 2) is the velocity in channel i. The Coulomb parameter in the second 
transfer channel is i] = 2M2{Z — l)/k'>/f> The coefficients K^" are obtained from the 

numerical solution of the Faddeev-Hahn-type equations. The cross sections are given by 

^ 47r(2L + l) {K-i)^ 

A;(i)2 {D - l)^ + {K^l + K^^)^ ^ ' 

where D = K^^K^i - K^2K^l When k^^^ 0, cx^u^^i ~ ^/k?. The transfer rates are 
defined by 

XXs^nl = <->nlVlNo (17) 

where vi is the relative velocity of the incident fragments and A^o the liquid-hydrogen density 
chosen here as 4.25 x 10^^ cm~^. We note that X^'^{k^^^ — 0) ~ const. In our model the total 
muon transfer rate is 

■^tot = X! -^Is^nl- (18) 
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III. NUMERICAL RESULTS 



We employ muonic atomic unit: distances are measured in units of a^, where is the 
radius of muonic hydrogen atom. The integro-differential equations are solved by discretizing 



them into a linear system of equations as in Refs. ||2T|J2^. As we are concerned with the 
low-energy limit, only the total angular momentum L = is taken into account. Even at 
zero incident energy, the transfer channels are open and their wave functions are rapidly 
oscillating Coulomb waves. A large number of discretization points is needed for achieving 
convergence. 

First Eqs. (|10D and (|l5) are solved setting the polarization potential to zero. The 
numerical convergence for the system of equations including the elastic and a transfer channel 
(at a time) is obtained after moderate effort. Finally, the total transfer cross section is 
calculated by adding the results of different two-channel contributions. In this case we need 
up to 700 discretization points per channel adequately distributed between and 70a^. 
This relatively large number of points is necessary to properly deal with the large final-state 
Coulomb interaction which could be as large as lle^/p in the case of argon. Some 20 to 
30 discretization points per channel would be enough for uncharged fragments in the initial 
and final state. 

Next the calculations are repeated in the presence of polarization potential. It is more 
difficult to obtain convergence with the polarization potential (|^ which is taken to be zero 
at small distances below the cut off A. In this case to get numerical convergence we have to 
integrate to very large distances — up to 300a^. We need up to 2000 discretization points 
per channel to obtain convergence. More points are employed near the origin and less at 
large distances. For example, near the origin we take up to 60 equally spaced points per 
unit length interval a^; in the intermediate region (p = 15 — 30a^) we take up to 8 equally 
spaced points per unit length interval, and in the asymptotic region (p = 30 — 300a^) we 
take up to 6 equally spaced points per unit length interval. 

The short-range potential of the present problem extends to about TZ = 20a^. It is 
customary to take the cut off A of the polarization potential much larger than the range TZ 
of the short-range potential. For A < TZ, the polarization potential becomes much larger 
than the short-range potential and for A = it diverges. The results for the cross section 
and the transfer rates can be very sensitive to the value of the cut off for A <TZ. However, 
the polarization potential should be effective in the intermediate region: oo > A > 7?.. We 
find that for oo > A > 7?-, the results for the transfer rates are very weakly dependent on A. 
We considered the polarization potential in this asymptotic region pi > A ~ 120a^. For a 
variation of A in this region from 120a^ to about 160a^, we find the transfer cross sections 
to be reasonably constant and the reported transfer rates of this study are the averages of 
these cross sections. If we increase A beyond 160a^ the effect of the polarization potential 
gradually decreases and transfer rates gradually tend towards the rates calculated without 
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the polarization potential. 

We present partial muon-transfer rates Ajijj and total transfer rates X\l^. We calculate the 
low-energy muon-transfer rates from {pfi)is and {dfi)is to Ne^°+, S^^"*", and Ar^^"*". We find 
in our calculation that in all cases the transfer takes place predominantly to the angular 
momenta states / = 0, 1 of the n = 6 orbital of (Ne^)^"*", n = 9 orbital of (S^)^^^, and n = 10 
orbital of (Ar^)^^+. The transition to other states of (X^)^^^^)+ is negligibly small. Hence in 
this work we only present muon-transfer rates for the / = 0, 1 states of the above orbitals of 
neon, sulphur, and argon. The low-energy partial rates A*[^^^„; /10^° s~^ and total rates A^q^ 
/10^° s~^ with and without the polarization potential are presented in Tables I to V together 
with the results of experimental works. In all cases the partial transfer rates without the 
polarization potential saturates to a reasonably constant value for E < 0.01 eV. A similar 
behavior is also observed in the presence of the polarization potential. In all cases presented, 
the rate of transition to the s state of the muonic atomic orbital is larger than that to the p 
state both in the presence and absence of polarization potential. For example, in the case of 
Ne^°^, the muon-transfer rate to the 6s state of (Ne^)^"*" is larger than that to the 6p state, 
and so on. 

First we consider the results in Table I for muon transfer from muonic protium {pfi)is 
to Ar^^"*". In this case we find that a constant value for the transfer rate is achieved for 
E < 0.04 eV and we present results up to this energy. The present total transfer rate of 
(12.9 ± 0.4) X 10^° s~^ is in good agreement with experiments of Refs. ||9Hri|. However, 



it disagrees strongly with experiments of Refs. |T^JT3[| . In the case of muon transfer from 
{dfi)is to Ar^^+, we find in Table II that the total transfer rate has decreased in this case 
compared to that in Table I, in agreement with the observed experimental trend For the 
rates reported in Table II a constant value is obtained for E < 0.04 eV. The present total 
transfer rate of (5.3 ± 0.4) x 10^° s^^ for {dfi)is to argon is in reasonable agreement with the 
experimental rates of Refs. @|14 . 



In the cases of muon-transfer rates from hydrogen isotopes to muonic sulphur and neon 
reported below in Tables HI, IV, and V, a constant value for the rate was not obtained till 
E = 0.04 eV, as in the case of argon above. Hence in these cases we report the results of 
our calculation for energies 0.01 and 0.04 eV with the experimental energy lying in between. 
In Table HI the results for muon transfer from muonic protium {pfJ^)is to S^^^ are shown. 
The present total transfer rate of (12.0 ± 0.4) x lO^'' s~^ is in reasonable agreement with 
the experimental rate 8.9 x 10^° s^^ of Ref. [0. In the case of muon transfer from muonic 
deuterium {dn)is to S^^^, we find from Table IV that the total transfer rate of (12.7±0.4) x 
10^° s~^ is also in good agreement with the experimental rate 11.0 x 10^° s~^ of Ref. |]T6[. The 
transfer rate has increased in this case compared to that in Table HI for muon transfer from 
protium in agreement with experimental observation [|16|. Finally, In Table V we exhibit 
muon transfer from muonic deuterium {dfi)is to Ne^°^. The present total transfer rate of 
(8.4 ±0.2) X 10^° s^^ at 0.01 eV is in good agreement with experimental rate 10.1 x 10^° s~^ 



9 



of Ref. m 



From Tables I to V we find that in most cases the transition rates calculated in the present 
model without the polarization potential is already in reasonable agreement with experiment. 
This suggests that the present model based on the close-coupling approximation to coupled 
Faddeev-Hahn-equations is very suitable for muon transfer reactions with neon, sulphur, and 
argon. A similar conclusion can be made from our previous study on muon transfer reactions 
with hydrogen [^, helium, lithium [^, carbon, and oxygen [^|. Although, the effect of 
polarization in these cases is expected to be large on the observables of the elastic-channel 
observables at low energies, this effect is not found to be so pronounced on the muon-transfer 
rates calculated in this work. 



IV. CONCLUSION 

We have calculated muon-transfer rates from muonic protium and deuterium to bare 
nuclei neon, sulphur, and argon employing a full three-body quantum-mechanical description 
of rearrangement scattering by solving the Faddeev-Hahn-type equations |]TP|] using close- 



coupling approximation. The exact inclusion of the final-state Coulomb interaction in the 
rearrangement channel (without approximation) guarantees the correct asymptotic form 
of the wave function in this channel. We also included a polarization potential at large 
distances in the initial channel beyond a cut off distance A. In these problems the short- 
range part of the interaction extends to about 20 muonic atomic units (25a^). By varying 
A in the intermediate range between 120a^ to 160a^ we find that the transfer rates are 
almost independent of the cut off A. This makes the reported transfer rates of this work 
quasi-independent of cut off. It is shown that in the present approach, the application of 
a close-coupling-type ansatz leads to satisfactory results for direct muon-transfer reactions 
from muonic hydrogen to bare nuclei neon, sulphur, and argon. The present rates are in 
good agreement with experiments p|- [lI| , p!4| -[T^ . Our calculation also reproduces the observed 
experimental trend of transfer rates to sulphur and argon from protium and deuterium 0. 
The transfer rate increases with the decrease of mass of the hydrogen isotope in the case of 
argon; the reverse is true for sulphur. 



It is interesting to note from the works of Refs. p(]|-[22[| and the present investigation that 
a low-order approximation to the Faddeev-Hahn-type equations as considered here produces 
very good muon-transfer rates from hydrogen isotopes to bare nuclei at low energies. In 
previous studies it was noted that low-order approximation to these equations also produces 
very accurate scattering observables at low energies for electron-hydrogen [^] and positron- 



hydrogen systems in agreement with the variational results. It is well-known that similar 
low-order approximation to the Schrodinger equation leads usually to poor results at low 
energies [p],^,|3l| . Hence low-order approximation to the Faddeev-Hahn-type equations is 



a very attractive alternative for studying low-energy Coulomb charge-transfer reactions in 
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general. 
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TABLES 



Table I. Low energy partial A^3_^„;/10^*^s ^ and total X^l^^/lO^^s ^ muon transfer rates reduced 



to liquid-hydrogen density A^^o = 4.25 x 10^^ cm ^ from muonic protium {pfi)is to hydrogen- 
excited state of muonic ar^ 
/IQio s- 



like excited state of muonic argon {ATf^)UJiQ together with some experimental results A*'' 



Energy 

E (eV) (n/) 



t^pol(Pl 



With polarization 



^tot 



Xtr 



Xtr 



Experiment 
Xtr 



0.04 10s 5.0 ±0.2 
lOp 3.9 ± 0.2 



8.1 ±0.2 
.9 ±0.4 4.8 ±0.2 



16.3 



12.0 Ig 14.6 ll 



12.9 ±0.4 



35 H 98 |13| 



0.06 10s 5.0 ±0.2 
lOp 3.9 ±0.1 



7.9 ±0.2 
.9 ±0.3 4.6 ±0.1 12.5 ±0.3 



0.1 10s 4.8 ±0.1 

lOp 3.8 ±0.1 



6.8 ±0.2 
.6 ±0.3 4.1 ±0.1 10.9 ±0.3 



Table 11. Same as in Table 1 from muonic deuterium {dfi)is to muonic argon (Ar^)^ 



17H 

mm=io 



Energy 
E (eV) 



(nl) 



f^pol(Pl 



With polarization 



Xtr 



xtr 

^tot 



Xtr 



Xtr 
^tot 



Experiment 
Xtr 



0.04 



10s 



1.3 ±0.1 
0.9 ±0.1 



2.2 ±0.2 



3.4 ±0.2 
1.9 ±0.2 



5.3 ±0.4 8.6 @ 9.4 |4 1 



0.06 



10s 
lOp 



1.2 ±0.1 
0.9 ±0.1 



2.1 ±0.2 



3.2 ±0.2 
1.8 ±0.2 



5.0 ±0.4 



0.1 



10s 
lOp 



1.1 ±0.1 
0.8 ±0.1 



1.9 ±0.2 



2.8 ±0.1 
1.4 ±0.1 



3.2 ±0.2 



Table 111. Same as in Table 1 from muonic protium {pfi)is to muonic sulphur (S^)„J9 



Energy 




f^pol(Pl) 


= 


With polarization 


Experiment 


E (eV) 


(nl) 






Xtr 




Xtr 


0.01 


9s 


6.5 




8.2 ±0.2 








9p 


3.1 


9.6 


3.8 ±0.2 


12.0 ±0.4 


8.9 [H 


0.04 


9s 


6.9 




8.4 ±0.2 








9p 


3.4 


10.3 


4.1 ±0.2 


12.5 ±0.4 
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Table IV. Same as in Table I from muonic deuterium (c?/i)is to muonic sulphur (S^)„J^ 



Energy 
E (eV) 


(nl) 


f^pol(Pl) 


= 
Xtr 


With polarization 

\tr \tr 
^Is^ni '^tot 


Experiment 


0.01 


9s 
9p 


6.8 

4.0 


10.8 


7.9 ±0.2 
4.8 ±0.2 


12.7 ±0.4 


11.0 


1 




0.04 


9s 

9p 


8.7 
4.4 


13.1 


9.7 ±0.2 
4.9 ±0.2 


14.6 ±0.4 








Table V. S 


ame as 


in Table I from muonic 


deuterium {dfj.)is to muonic neon (Ne^)^lg 


Energy 
E (eV) 


(nl) 


t^pol(Pl) 


= 
Xtr 


With polarization 

\tr \tr 


Experiment 
Xtr 


0.01 


6s 
6p 


6.9 
< 0.1 


7.0 


8.2 ±0.2 
< 0.2 


8.4 ±0.2 


10.1 1 


15 


1 


0.04 


6s 
6p 


4.7 
< 0.05 


4.75 


6.2 ±0.2 
< 0.1 


6.3 ±0.2 
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